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Cosmos  236  (1968-70A)  was  launched  on  27  August  1968  into  a  near-circular 
orbit  of  inclination  56f  and  is  expected  to  decay  during  late  1989.  The  orbit 
has  been  determined  from  observations  for  77  epochs  between  July  1983  and 
October  1984  over  the  time  interval  when  the  orbit  was  expected  to  be  signifi¬ 
cantly  influenced  by  the  effects  of  15th-order  resonance  with  the  Earth's  gravi¬ 
tational  field:  exact  resonance  occurred  on  13  March  1984.  The  observations 
numbered  over  4700,  including  284  from  the  Hewitt  cameras -of  the  University  of 
Aston 'which  are  sited  at  Herstmonceux  in  Eny^and  and  at  Siding  Spring  in 
Australia.  The  orbital  accuracy  achieved  was  fairly  consistent  throughout,  with 
the  standard  deviation  in  orbital  inclination  and  eccentricity  corresponding  on 
average  to  positional  accuracies  of  85  m  and  65  ra  respectively. 

Analysis  of  the  changes  in  inclination  and  in  eccentricity  at  resonance 
has  given  values  of  three  pairs  of  lumped  harmonics  of  order  15  and  three  pairs 
of  order  30,  one  pair  of  each  from  inclination  and  two  from  eccentricity.  The 
values  from  inclination  had  standard  deviations  equivalent  to  accuracies  in 
geoid  height  of  0.6  cm  and  2.0  cm  for  orders  15  and  30  respectively  while  the 
equivalent  accuracies  for  the  values  from  eccentricity  were  1.6  cm  and  6.0  cm. 

)l 's  jtff-.  j  3  -V  ") 

i 

Departmental  Reference:  SS  9 


Copyright 

© 

Controller  HMSO  London 
1988 


UNLIMITED 


2 


LIST  OF  CONTENTS 

Page 


1  INTRODUCTION  3 

2  THE  OBSERVATIONS  AND  THE  ORBITS  3 

2.1  Sources  of  the  observations  3 

2.2  Observational  accuracy  3 

2.3  The  orbits  5 

3  RESONANCE  THEORY  6 

4  THE  PASSAGE  THROUGH  RESONANCE  FOR  COSMOS  236  7 

5  ANALYSIS  OF  INCLINATION  7 

6  ANALYSIS  OF  ECCENTRICITY  8 

7  INCLINATION  AND  ECCENTRICITY  FITTED  SIMULTANEOUSLY  9 

8  THE  EFFECT  OF  THE  DAY-TO-NIGHT  VARIATION  IN  AIR  DENSITY  10 

9  THE  ACCURACY  OF  THE  LUMPED  HARMONICS,  IN  TERMS  OF  GEOID  HEIGHT  11 

10  CONCLUSIONS  12 

Appendix  Values  of  Q  coefficients  for  1968-70A  15 

Table  2  17 

References  21 

Illustrations  Figures  1-5 


Report  documentation  page 


inside  back  cover 


058 


058 


3 


t  INTRODUCTION 

The  satellite  Cosmos  236,  designated  1968-70A,  was  launched  on 
27  August  1968  and  is  expected  to  decay  in  the  Earth's  atmosphere  during  the  last 
half  of  1989.  It  is  cylindrical  in  shape,  about  2  m  long  and  1  m  diameter,  with 
a  weight  of  about  850  kg  and  its  initial  orbital  parameters'  were:  inclination 
56.07°,  perigee  and  apogee  heights  588  and  630  km  respectively,  and  nodal  period 
96.83  minutes. 

The  orbit  of  Cosmos  236  contracted  slowly  under  the  influence  of  air  drag, 
and  in  March  1984  it  passed  through  15th-order  resonance,  which  occurs  when  a 
satellite's  track  over  the  Earth  repeats  after  15  revolutions,  is  when  the 
satellite  makes  15  revolutions  while  the  Earth  spins  once  relative  to  the  orbital 
plane.  If  the  passage  through  resonance  is  slow  enough,  the  effects  of  15th-crder 
harmonics  in  the  geopotential  can  build  up  and  result  in  an  appreciable  pertur¬ 
bation  to  some  of  the  orbital  elements:  this  variation  can  be  analysed  to  derive 
values  for  lumped  geopotential  harmonics  of  order  15.  The  orbit  of  Cosmos  23b 

has  been  determined  between  July  1983  and  October  1984  from  radar  and  optical 
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observations  with  the  aid  of  the  RAF.  orbit  refinement  program  PRCP6,  and  the 
changes  in  inclination  and  eccentricity  have  been  analysed  to  give  six  values  of 
lumped  15th-ordcr  harmonics  and  six  of  30th-order. 

2  THE  OBSERVATIONS  AND  THE  ORBITS 
2  .  1  Sources  of  the  observations 

The  orbit  of  1968-70A  has  been  determined  at  77  epochs  between  4  July  1983 
and  30  October  1984  from  4744  observations,  not  including  those  rejected  in  the 
orbit  determinations. 

The  observations  came  from  three  different  sources,  the  most  accurate  being 
those  from  the  University  of  Aston's  Hewitt  cameras  at  the  Royal  Greenwich 
Observatory,  Herstmonceux,  and  at  Siding  Spring  in  Australia;  284  of  these 
observations  were  used  in  33  of  the  77  orbits.  The  second  group  consisted  of 
282  visual  observations  made  by  volunteer  observers  reporting  to  the  Earth 
Satellite  Research  Unit  at  the  University  of  Aston.  The  third  and  largest  group 
were  4178  radar  observations  made  by  the  US  Navy  Navspasur  system,  kindly  supplied 
by  the  US  Naval  Research  Laboratory. 

2 . 2  Observational  accuracy 

The  rms  residuals  of  the  observations  have  been  calculated  using  the  RAE 
3 

computer  program  ORES  and  have  been  sent  to  the  observers.  Table  1  gives  the 
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residuals  for  selected  observing  stations  with  at  least  six  observations  accepted 
in  the  final  orbit  determinations.  The  US  Navy  observations  are  geocentric,  and 
if  they  were  given  in  the  same  form  as  the  topocentric  observations,  their  angular 
rms  residuals  would  increase  by  a  factor  of  between  5  and  10.  In  calculating  the 
rms  residuals  for  the  visual  observers,  observations  with  residuals  greater  than 
twice  the  rms  have  been  omitted,  the  numbers  used  being  shown  in  brackets.  This 
gives  a  truer  impression  of  the  normal  accuracy  of  the  observer,  as  it  eliminates 
observations  marred  by  poor  seeing  conditions  and  possible  deficiencies  in  orbital 
fitting. 


Table  1 


Residuals  for  selected  stations 


Station 

Number  of 
accepted 
observations 

rms  residuals 

Range 

km 

Minutes  of  arc 

RA 

Dec 

Total 

1 

US  Navy 

565 

2.2 

1  .7 

2.8 

2 

US  Navy 

338 

2.9 

3.0 

4.2 

3 

US  Navy 

395 

3.1 

2.4 

3.9 

4 

US  Navy 

533 

3.1 

2.7 

4.2 

5 

US  Navy 

369 

2.6 

2.2 

3.4 

6 

US  Navy 

554 

2.3 

1  .9 

3.0 

29 

US  Navy 

1424 

0.6 

0.3* 

0.2* 

414 

Capetown 

12(11) 

1  .7 

1.7 

2.4 

2122 

Malvern  5 

11 

2.0 

1  .6 

2.6 

2265 

Farnham 

6(5) 

3.2 

2.5 

4.0 

2414 

Bournemouth 

112(104) 

4.0 

4.2 

5.8 

2418 

Sunningdale 

15 

4.1 

4.5 

6.1 

2420 

Willowbrae 

54(51) 

7.9 

3.9 

8.8 

2437 

Warrington 

12 

5.8 

6.3 

8.6 

2539 

Dymchurch 

15(14) 

2.0 

1  .6 

2.6 

2657 

Bridgwater 

10 

1  .6 

2.4 

2.9 

2659 

Herstmonceux  3  (Hewitt  camera) 

266(250) 

0.11 

0.06 

0.13 

4156 

Apeldoorn 

14(12) 

3.1 

2.6 

4.1 

9652 

Siding  Spring  (Hewitt  camera) 

18 

0.04 

0.05 

*  Geocentric 


The  rms  residuals  of  the  Hewitt  cameras  are  8  seconds  of  arc  from  250 
observations  by  the  Herstmonceux  camera  and  3  seconds  of  arc  from  18  observations 
by  the  Siding  Spring  camera,  equivalent  to  about  30  and  10  m  in  position  respec¬ 
tively.  The  accuracy  in  plate  reading  is  usually  estimated  as  about  2  seconds  of 
arc  and  the  residuals  for  Siding  Spring  confirm  this  estimate.  The  residual  is 
of  course  a  combination  of  (a)  observational  error  and  (b)  errors  in  the  orbital 
model  and  the  orbit  determination  process.  A  number  of  the  orbits  utilized  more 
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than  one  plate  from  the  Herstmonceux  camera,  and  the  PROP  orbital  model,  which 
ignores  lunisolar  perturbations,  is  not  accurate  enough  to  do  justice  to  the 
multiple  plates  in  one  orbit  determination:  lunisolar  perturbations  can  build 
up  to  the  order  of  50-100  m  in  a  week,  much  greater  than  the  5-10  m  capability 
of  the  cameras.  The  Herstmonceux  camera  residuals  are  therefore  as  good  as  could 
be  expected  with  PROP. 

2 . 3  The  orbits 

The  orbits  have  been  determined  with  the  aid  of  the  RAE  orbit  refinement 
program  PROP  in  the  PR0P6  version,  and  the  orbital  elements  at  the  77  epochs, 
with  their  standard  deviations,  are  given  in  Table  2  (page  17).  The  epoch  for 
each  orbit  is  at  00  hours  on  the  day  indicated,  and  the  PROP  program  fits  the 
mean  anomaly  M  by  a  polynomial  of  the  form 

M  =  M.  +  M  t  +  M„t2  +  M,t3  +  M. tA  +  M,t5  (1) 

0  12  3  4  5 

where  t  is  the  time  measured  from  epoch.  Up  to  six  M-coef f icients  may  be  used, 
depending  on  the  severity  of  the  drag.  For  1968-70A,  which  at  resonance  was  in 
a  near-circular  orbit  at  a  height  of  about  500  km,  Mq  -  were  sufficient  for 
59  of  the  orbits,  and  18  orbits  required  . 

All  of  the  orbits  fitted  the  observations  satisfactorily,  with  the  value 
of  t  ,  the  parameter  which  indicates  the  measure  of  fit,  varying  from  0.26  to 
0.90  with  an  average  value  of  0.53.  The  average  number  of  observations  in  an 
orbit  determination  was  62,  spread  over  a  time  interval  averaging  5.9  days. 

The  average  standard  deviation  in  eccentricity  for  all  the  77  orbits  is 
0.000009,  equivalent  to  an  error  in  perigee  distance  of  65  m,  and  this  is  much 
the  same  as  the  average  standard  deviation  for  the  33  orbits  containing  Hewitt 
camera  observations.  As  the  eccentricity  is  very  small,  it  is  useful  to  plot 
the  values  in  polar  form,  as  shown  in  Fig  1 ,  where  the  values  are  linked  with 
a  continuous  line  as  a  guide  to  the  eye.  In  the  absence  of  drag  and  resonance, 
the  locus  of  the  points  should  be  close  to  a  circle;  when  drag  acts,  the  circle 

It  .  . 

is  converted  to  a  contracting  spiral  .  lu  Fig  1,  however,  the  initial  circuit 
(orbits  1-29)  is  followed  by  a  second  (orbits  30-57)  which  spirals  outwards, 
presumably  as  a  result  of  the  resonance.  (Thus  it  might  be  guessed,  even  from 
the  raw  values  of  Fig  1 ,  that  the  resonance  begins  to  act  strongly  at  about 
orbit  29.)  The  third  circuit  closely  follows  the  first  between  orbits  61  and  71. 
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The  mean  standard  deviation  in  inclination  for  the  77  orbits  is  0.0007°, 
corresponding  to  an  error  of  about  85  m  in  cross-track  distance;  for  the  33  orbits 
with  Hewitt  camera  observations  the  accuracy  is  better,  the  average  standard 
deviation  being  0.0005°. 

3  RESONANCE  THEORY 

The  theory  has  often  been  given  (eg  Ref  5)  and  will  only  be  outlined  here. 

In  brief,  the  rate  of  change  of  inclination  i  caused  by  a  relevant  pair  of 

geopotential  coefficients  C„  and  S„  near  resonance  may  be  written 

Jim  Jim 


di 

dt 


n\1  -el 

f-1 

sin  i 

F„  G.  (k 
Jimp  Hpq 


cos  i  - 


'  j§fm)  eXP{j(Y<t  - 


.  (2) 

where  F  and  G  are  functions  of  inclination  and  eccentricity  defined  in  Ref  5, 
R  is  the  Earth's  equatorial  radius,  and  the  resonance  angle  4>  for  15:1 
resonance  is  given  by 


$  =  u  +  M  +  15(Q  -  v) 


(3) 


where  v  is  the  sidereal  angle.  The  indices  y  and  q  are  integers,  and  in 

practice  the  most  important  terms  are  those  with  y  =  1  ,  though  those  with 

y  =  2  may  also  be  needed.  For  inclination,  the  q  =  0  terms  are  the  most 

important,  but  the  eccentricity  is  affected  most  by  the  terms  with  q  =  ±1  . 

For  15:1  resonance  the  equations  linking  m,  y,  q,  k,  p  and  1.  are:  m  =  15y  ; 

k  =  y-  q;  2p  =  £-  k.  The  values  of  l  must  be  such  that  SL  ^  m  and  (Jl  -  k) 

is  even.  The  successive  C„  and  S„  coefficients  that  arise  may  be  grouped 

Jim  Jim 

into  a  'lumped  harmonic'. 


Jl 


(4) 


where  Jl  increases  in  steps  of  2  from  its  minimum  value 

serve 
_q,k 


are 


0  -  "‘S’ 

constant  coefficients.  Thus  the  observed  change  in  i  at  resonance  gives  values 

and  S 


and  the 

ve 

appropriate  to  a  particular  inclination, 


_q,k 

of  the  lumped  harmonics  C  _  _ 

mm 

here  56.08  .  When  values  are  available  for  many  different  inclinations,  the 
individual  harmonics  can  be  determined^. 


The  rate  of  change  of  eccentricity  produced  by  the  (l,m)  harmonics  is 

,  5 

given  by  : 
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de 

dt 


£pq 


q  -  j (k  +  3q)e2 
e 


.  L-m+ 1 
J 


jS2m)  exp  j(Y$  -  qo>) 


. (5) 

Again  the  C.  and  S.  may  be  grouped  into  appropriate  lumped  harmonics. 

stm  V-m 

4  THE  PASSAGE  THROUGH  RESONANCE  FOR  COSMOS  236 

Exact  resonance ,  i  =  C  ,  occurred  on  13  March  1984,  and  the  variation  of 
$  and  $  with  time  is  shown  in  Fig  2.  It  can  be  seen  that  $  increases  fairly 
steadily  between  -6.8  degrees/day  initially  and  +5.3  degrees/day  at  the  end,  so 
the  resonance  is  well  balanced  and  does  not  suffer  from  the  considerable  changes 
in  drag  which  affect  many  resonance  analyses.  This  good  balance  is  due  to  the 
fortunate  chance  that  resonance  occurred  at  a  time  close  to  the  minimum  of  the 
solar  cycle,  when  air  density  was  at  its  lowest  and  fairly  steady. 

5  ANALYSIS  OF  INCLINATION 

Two  important  perturbations  need  to  be  subtracted  from  the  raw  values  of 
inclination  in  Table  2.  The  first  is  that  due  to  the  combined  effect  of  lunisolar 
and  zonal  harmonic  perturbations,  which  were  calculated  using  the  RAE  computer 
program  PROD7,  with  a  1-day  integration  interval  and  restarts  every  20  days  or 
less.  The  second  is  that  due  to  the  J,  ,  tesseral  harmonic,  which  is  printed 
on  each  orbit  determination  output  by  the  PROP  program.  These  two  perturbations, 
which  had  maximum  numerical  values  of  0.0024°  and  0.0017°  respectively,  were 
subtracted  from  the  raw  values.  Perturbations  due  to  tides  and  solar  radiation 
pressure  were  considered  too  small  to  be  worth  evaluating.  The  resulting  values 
of  inclination,  with  standard  deviations,  are  shown  in  Fig  3:  the  main  change  due 
to  resonance  and  the  subsidiary  oscillations  are  well  displayed. 

The  values  were  then  fitted  with  an  integrated  form  of  the  theoretical 

g 

equation  (2)  with  the  aid  of  the  THROE  computer  program  ,  for  various  pairs  of 

9 

(y,q)  >  assuming  an  atmospheric  rotation  rate  of  0.9  rev/day;  a  minimum  standard 
deviation  of  0.0005°  was  set,  in  view  of  the  neglect  of  earth  and  ocean  tides. 

The  THROE  program  adjusts  the  values  of  i  for  the  effects  of  atmospheric 
rotation  and  lunisolar  precession  of  the  Earth's  axis. 

As  a  result  of  previous  analyses  of  inclination1^  for  near-circular  orbits 
at  15th-order  resonance,  it  was  expected  that  the  most  important  pairs  of  values 
of  (y,q)  would  be  (y,q)  =  (1,0)  and  (2,0)  with  the  possible  additional  pairs 
(y,q)  “  (1,1)  and  (1,-1).  Here  it  was  found  that  there  was  no  advantage  in 


00 
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using  the  additional  pairs  because  the  fitting  was  not  substantially  improved  and 
the  lumped  values  derived  were  indeterminate.  For  similar  reasons  the  addition 
of  (Y,q)  =  (3,0)  terms  was  unavailing.  The  fitting  was  therefore  made  with 
(y,q)  =  (1,0)  and  (2,0)  only:  the  nominal  standard  deviations  were  doubled  on 
eight  orbits*  to  avoid  residuals  of  more  than  twice  the  measure  of  fit  e  ,  and 
one  standard  deviation  (for  MJD  45872)  had  to  be  quadrupled,  for  the  same  reason. 


The  THROE  fitting  of  the  theoretical  curve  to  the  observed  values  is  shown 
in  Fig  4;  the  orbits  mentioned  above  are  shown  with  their  relaxed  standard 
deviations.  The  measure  of  fit  e  had  the  value  1.49  and  the  values  of  the 
lumped  harmonics  obtained  were  as  follows: 


9-  » 

10*C15  -  -213.5  ±  5.4 

q.0.2 

10  C3Q  =  -34  ±  149 

The  fitting  in  Fig  4  is  very  good,  and 
The  value  of  e  substantially  greater 
deviations  are  slightly  overoptimistic 
imperfections . 


Q_V  ,  i 

107S15  =  -91  .  i  ±  4.2 

q_0,2 

10  S30  =  -624  ±  106  •  (6) 

the  values  derived  should  be  reliable, 
than  1  implies  either  that  the  standard 
or  that  the  modelling  has  some  small 


6  ANALYSIS  OF  ECCENTRICITY 

The  lunisolar  perturbations  to  eccentricity  e  have  been  obtained  using 
the  PROD  program,  and  have  been  removed.  The  air  drag  effects  have  been  removed 
within  THROE,  assuming  a  constant  scale  height  H  of  60  km,  appropriate  to  a 
height  of  490  km,  and  taking  mean  values  of  between  successive  orbits''. 

The  zonal  harmonic  perturbations  were  also  removed  within  THROE  and  it  was  found 
that  a  small  adjustment  to  the  odd  zonal  harmonic  oscillation,  expressed  as  a 
change  in  10  J3  from  -2.53  to  -2.60,  was  beneficial  in  improving  the  fitting. 

The  values  of  e  thus  obtained  are  shown  in  Fig  5,  which  also  gives  the 
final  THROE  fitting,  with  (y,q)  =  (1,1),  (1,-1),  (2,1)  and  (2,-1).  It  is 
immediately  apparent  in  Fig  5  that  the  THROE  fitting  is  deficient  between 
MJD  45782  and  45818,  orbits  44-50,  and  Fig  1  shows  that  this  is  just  the  region 
where  e  becomes  extremely  small  and  iL  varies  greatly  (w  decreases  by  115° 
in  5  days  between  orbits  50  and  51).  THROE  is  not  fully  adapted  for  avoiding 
problems  when  e  is  very  small,  because  it  is  formulated  in  terms  of  e  and  ou  , 


*  At  MJD  45519,  45523,  45531,  45647,  45654,  45751,  45805,  45978. 
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rather  than  e  cos  w  and  e  sin  u  ,  as  is  necessary  for  very  small  e  .  Thus 
the  reliability  of  the  results  from  this  fitting  is  open  to  question,  although 
it  may  be  that  the  high  value  of  the  measure  of  fit  will  in  itself  provide 
sufficient  allowance  for  the  ill-fitting  region.  In  the  fitting  the  seven  values 
between  MJD  45782  and  45818  were  all  relaxed,  the  first,  sixth  and  seventh  by  a 
factor  of  2  and  the  others  by  a  factor  of  4,  to  avoid  residuals  >2e  .  Apart 
from  this  'problem  region',  it  was  only  necessary  to  relax  three  values  by  a 
factor  of  2  (MJD  45536,  45597  and  45654)  and  one  by  a  factor  of  4  (MJD  45527'. 
With  these  standard  deviations,  which  are  indicated  in  Fig  5,  the  measure  oi  fit 
e  was  2.63  and  the  values  of  the  lumped  harmonics  obtained  were: 


9  '.0 

10  C,5 

CM 

+! 

CM 

II 

9_1  ,0 
,0  s15 

=  39.3  ±  7.6 

9--’>2 

.0  c,5 

=  74.4  ±  5.8 

-1  2 

109S  ’ 

1 } 

=  37.5  ±  4.4 

o  J.1 
to  c30 

=  -343  ±  121 

9J.> 

10  s30 

-  -116+234 

9--'  ■  J 

to  c30 

-  553  +  77 

—  o 

1  CO 

I'C 

O' 

o 

-  509  ±  117 

{  7 

In  several  recent  analyses  of  eccentricity  at  l5th-order  resonance  the 
effect  of  the  day-to-night  variation  in  density  has  been  taken  into  account. 

Here  this  effect  is  quite  small,  as  indicated  in  section  8,  and  has  not  been 
taken  into  account. 

7  INCLINATION  AND  ECCENTRICITY  FITTED  SIMULTANEOUSLY 

The  changes  in  inclination  ; ^d  eccentricity  have  been  fitted  simultaneously 
with  the  aid  of  R.H.  Gooding's  program  SIMRES.  It  was  not  expected  that  better 
values  would  emerge,  because  the  sets  of  lumped  harmonics  obtained  from  i  are 
completely  different  from  those  obtained  from  e  ;  but  the  simultaneous  fitting 
seemed  worth  trying. 

As  usual,  the  inclination  and  eccentricity  were  weighted  in  accordance  with 
the  t  values  of  the  contributing  THROE  runs,  and  the  values  obtained  for  the 
six  pairs  of  lumped  harmonics  are: 


i 
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9-°’ 1 
.0  C,5 

=  -211.115.2 

Q_0,1 

1°  s15 

=  -87.3  ±  4.0 

q_0-2 

io  c3Q 

=  65  ±  143 

q-0. 2 

10  S30 

=  -624  ±  104 

oj,0 
IOC, 3 

=  19.3  1  10.4 

q  1,0 

10  S  _ 

1  D 

=  39.3  1  7.7 

9_-’ >2 
,°c,5 

=  74.8  ±  6.0 

9__1 >2 
.0  S,5 

=  37.0  ±  4.5 

c  1,1 

,o  c30 

=  -347  ±  123 

o_1,' 

10  S30 

=  -114  ±  238 

q_-1,3 
i°  c3Q 

=  559  1  78 

9_” 1 , 3 

,0S30 

»  507+118 

It  will  be  seen  that  the  last  four  pairs  of  values  are  virtually  the  sane 

as  in  equations  (7)  derived  from  eccentricity  alone,  though  the  standard  deviations 

are  marginally  greater:  this  is  to  be  expected,  because  the  first  two  pairs  of 

values  have  very  little  influence  on  the  eccentricity.  The  first  two  pairs  of 

values  are,  however,  appreciably  different  from  those  obtained  from  inclination 

.0,1 

alone,  equations  (6),  the  largest  change  (for  )  being  0.9  sd.  Again  this 

is  to  be  expected,  because  the  other  four  pairs,  the  values  of  which  are  in 

effect  'imposed'  by  e  ,  do  have  a  slight  influence  on  i  .  If  the  values 

obtained  from  e  were  fully  reliable,  the  SIMRES  solution  would  be  preferred; 

as  the  e-analysis  is  slightly  questionable,  however,  the  separate  solutions, 

equations  (6)  and  (7),  have  to  be  recommended.  (It  is  therefore  a  little  ironical 

6  -0,1 

that  tn  the  solutions  for  individual  coefficients  ,  the  SIMRES  values  of  S... 

_0,2  _  .0,1  15 

and  C^q  fit  better,  though  fits  worse.) 

8  THE  EFFECT  OF  THE  DAY-TO-NIGHT  VARIATION  IN  AIR  DENSITY 

The  day-to-night  variation  in  air  density  has  an  effect  on  eccentricity,  but 
has  been  neglected  and  needs  to  be  approximately  estimated.  When  e  <  0.003  , 
as  it  is  here,  the  decrease  in  e  per  revolution  due  to  drag  is  given  by 
equation  (11.11)  of  Ref  4  as 

Ae  =  -itiar,^  exp|(rQ  -  a)/H^  ^  z  +  F  cos  $  +  O^e,  -r  Z^)J  .  (9l 
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where  6  is  the  drag  parameter,  the  density  at  distance  rQ  when  $  =  90  , 

and  z  =  ae/H  =  !15e  here,  with  H  =  60  km  ,  so  that  z  is  generally  of  order 
0.3.  The  term  F  cos  .)>  in  equation  (9)  expresses  the  effect  of  the  day-to¬ 
night  variation  in  density:  here  the  appropriate  value  of  F  is  about  0.6  (see 
Table  7.1  of  Ref  4),  and  the  geocentric  angle  <J>  between  the  perigee  and  the 
centre  of  the  'diurnal  bulge'  goes  through  six  cycles.  As  a  first  approximation, 
therefore,  ignoring  the  F-term  in  (9)  is  equivalent  to  ignoring  a  six-cycle 
oscillation  having  an  average  amplitude  twice  as  large  as  the  spherical-atmosphere 
term  (because  F  ~  2z) .  The  maximum  effect  caused  by  the  day-to-night  variation 
is  2/r  times  twice  the  spherical-atmosphere  effect  over  a  half-cycle  of  ;  , 

which  is  about  one-twelfth  of  the  total  time  interval.  The  spherical-atmosphere 
correction  to  e  ,  which  has  been  evaluated  within  THROE,  is  almost  linear  and 
reaches  0.000143  at  the  final  orbit;  thus  it  is  about  0.000012  during  one-twelfth 
of  the  time  interval.  Therefore  the  maximum  change  in  e  caused  by  the  day-to- 
night  effect  is  on  average  0.000015  (on  applying  the  factor  4/r).  This  assumes 
that  <j>  ^  goes  through  a  full  cycle,  but  in  fact  the  average  minimum  of  i  in 
the  six  cycles  is  45*'  and  the  average  maximum  140°,  so  tne  average  amplitude  of 
the  oscillation  is  about  252  less  than  estimated  above,  namely  0.00001 1.  This 
is  less  than  the  average  sd  of  the  values,  with  the  relaxations  used  in  Fig  5: 
since  r.  -  2.63  ,  the  inclusion  of  the  effects  of  the  day-to-night  variation 
would  be  unlikely  to  have  a  significant  effect. 


However,  it  sr  happens  that  between  MJD  45766  and  45805,  at  the  start  of 
the  ill-fitting  region  in  Fig  5,  the  perigee  is  continually  within  90°  of  the 
centre  of  the  bulge,  and  the  day-to-night  effect  is  unusually  large,  increasing 
to  about  0.000030.  However,  the  points  between  MJD  45790  and  45814  demand  a 
correction  of  0.000150  to  bring  them  close  to  the  curve;  so  even  this  extra- 
large  correction  is  much  less  than  required,  though  it  would  slightly  improve 
the  fitting. 

It  can  be  concluded  that  the  fitting  in  Fig  5  would  not  be  much  improved 
by  taking  account  of  the  day-to-night  variation  in  air  density.  If,  in  future, 
a  better  analysis  can  be  made  through  improvements  in  the  THROE  program,  the 
day-to-night  correction  would  be  worth  making. 


9  THE  ACCURACY  OF  THE  LUMPED  HARMONICS,  IN  TERMS  OF  GEOID  HEIGHT 

The  standard  deviations  o  of  the  values  of  lumped  harmonics  in  equations 
(6)  and  (7)  can  be  approximately  interpreted  as  equivalent  to  accuracies  "  , 

say,  in  geoid  height.  The  linking  equation  is  0^  ~  Ro/Q  , 
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where  Q  =  I  £  (<%* 'I' rf }  * 

l 

on  the  assumption  that  the  magnitudes  of  the  individual  coefficients  fall  off  as 
2 

1/£  .  The  numerical  values  of  the  Q  coefficients  for  1968-70A  are  listed  in 

the  Appendix,  and  the  values  of  o^  for  the  lumped  harmonics  in  equations  (6) 
and  (7)  are  as  follows. 


Table  3 

for  15th-order  lumped  harmonics 


1 

0,1 

0,1 

1,0 

1,0 

.-1,2 

-1,2 

;  Harmonic 

C,5 

S,5 

C,5 

S 1 5 

C1 5 

S15 

1 

!  o  (cm) 

!  8 

0.7 

0.55 

1.8 

1.3 

2.0 

1  .5 

Table  4 

o  for  30th-order  lumped  harmonics 


0,2 

0,2 

.1,1  | 

.1,1 

.-',3 

-',3 

Harmonic 

C30 

S30 

So 

So 

So 

So 

n  (cm) 

g 

B 

1.7 

B 

B 

5.3 

8.1 

Thus  it  appears  that  the  15th-order  lumped  harmonics  determined  from  the 
inclination  have  accuracies  equivalent  to  about  0.6  cm  in  geoid  height,  while 
those  determined  from  eccentricity  have  an  average  accuracy  of  1 .6  cm  in  geoid 
height.  For  the  30th-order  lumped  harmonics  the  corresponding  accuracies  are 
2  cm  and  6  cm. 

10  CONCLUSIONS 

The  orbit  of  1968-70A  was  significantly  influenced  by  the  effects  of  15th- 
and  30th-order  harmonics  in  the  geopotential  during  all  the  15  months  over  which 
the  orbit  was  determined,  from  4  July  1983  to  30  October  1984.  The  orbit  deter¬ 
minations  were  at  77  epochs,  and  utilized  4744  observations  including  284  Hewitt 
camera  observations:  the  accuracy  achieved  was  fairly  uniform,  the  average 
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standard  deviations  in  inclination  and  eccentricity  being  about  85  m  and  65  m 
respectively.  This  is  the  most  accurate  and  most  extensive  orbit  to  be  analysed 
at  15th-order  resonance  at  an  inclination  near  56°. 

Analysis  of  the  variation  in  orbital  inclination  (Fig  4)  has  yielded  values 
of  lumped  harmonics  of  order  15  and  30,  given  in  equations  (6),  which  have 
standard  deviations  equivalent  to  accuracies  in  geoid  height  of  0.6  cm  and  2.0  cm 
respectively.  Analysis  of  the  variation  in  orbital  eccentricity  gave  values  of 
two  pairs  of  lumped  harmonics  of  order  15  and  two  of  order  30:  see  equations  (7). 
The  fitting  of  eccentricity  was  less  than  perfect  (see  Fig  5),  but  the  accuracies 
were  still  good,  being  equivalent  to  an  average  of  1.6  cm  in  geoid  height  for 
15th  order  and  6  cm  for  30th  order.  The  results  have  already  been  used^  to 
improve  the  determination  of  the  individual  harmonic  coefficients  of  order  15 
and  30. 


Appendix 


VALUES  OF  THE  Q  COEFFICIENTS  FOR  1968-70A 
Q  k 

The  Q  coefficients  Q?’  determine  the  numerical  dependence  of  the 
q  k  ^ 

lumped  harmonic  C  on  the  indi"idual  harmonic  coefficients  C„  ,  as  indi 

m  Jim 

cated  in  equation  (4).  For  m  =  15  and  (q,k)  =  (0,1)  ,  for  example, 

^ '  -  n  i-  n 


r  +  n0*1?  +  n0’1? 

15,15  Q17  C 1 7  , 1 5  yi9  19,15 


For  m  =  15  and  (q,k)  =  (1,0)  ,  the  first  term  has  £,  even,  and 

_ ^  -  i  i  r»_ 


?  +  n’-°F  +  n1 »°r 

16,15  ^18  L18,15  V20  20,15 


The  equations  are  similar  for  S  ,  substituting  S  for  C  .  The  values  of 
the  relevant  Q  coefficients  are  tabulated  below. 

Table  5 

Values  of  for  15th  order  (m  =  15) 


(<l,k)  =  (0,1) 


Q 


(q,k)  =  (1,0)  (q.k)  =  (-1,2) 


Append 


Table  6 

Values  of  for  30th  order  (m  =  30) 


,k)  =  (0,2) 

(q,k)  -  (1,1) 

_ 

X) 

V 

II 

1 

u> 

Q 

i 

Q 

l 

Q 

1  .000 

31 

1  .000 

31 

1  .000 

-8.756 

33 

-5.683 

33 

-4.362 

26.447 

35 

14.752 

35 

8.045 

-38.430 

37 

-19.658 

37 

-5.902 

21.820 

39 

9.670 

39 

-1  .795 

10.542 

41 

7.296 

41 

4.672 

-17.172 

43 

-9.471 

43 

0.671 

-4.024 

45 

-3.231 

45 

-3.392 

11.453 

47 

7.066 

47 

-0.838 
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inclination  after  removal  of  lunisolar,  zonal  harmonic  and  J2  2  perturbations 


Fig  4  Values  of  inclination  with  fitted  curve  for  (7, q)  =  (1,0),  (2,0) 
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1  7  \bstr.iet  , 

Cosmos  236  (1963-70A)  was  launched  on  27  August  1968  into  a  near-circular 
orbit  of  inclination  56°  and  is  expected  to  decay  during  late  1989.  The  orbit  has 
been  determined  from  observations  for  77  epochs  between  July  1983  and  October  1984 
over  the  time  interval  when  the  orbit  was  expected  to  be  significantly  influenced  by 
the  effects  of  S5th-order  resonance  with  the  Earth's  gravitational  field:  exact 
resonance  occurred  on  13  March  1984.  The  observations  numbered  over  4700,  including 
284  from  the  Hewitt  cameras  of  the  University  of  Aston  which  are  sited  at 
llurstmonccux  in  England  and  at  Siding  Spring  in  Australia.  The  orbital  accuracy 
achieved  was  fairly  consistent  throughout,  with  the  standard  deviation  in  orbital 
inclination  and  eccentricity -c.c ‘responding  on  average  to  positional  accuracies  of 
85  m  and  65  m  respectively. 

Analysis  of  the  changes  in  inclination  and  in  eccentricity  at  resonance  has 
given  values  of  three  pairs  of  lumped  harmonics  of  order  15  and  three  pairs  of 
order  30,  one  pair  of  each  from  inclination  and  two  from  eccentricity.  The  values 
from  inclination  had  standard  deviations  equivalent  to  accuracies  in  geoid  height 


>f  0.6  <  rn 


1.0  cm  for  orders  15  and  30  respectively  while  the  equivalent 


f  or  the  values  from  eccentricity  were.  1  .6' cm  and  6.0  cm. 
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